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In this paper, the performance of cognitive transmission under quahty of service (QoS) constraints and interference 

hmitations is studied. Cognitive secondary users are assumed to initially perform sensing over multiple frequency 

Q>^ ' bands (or equivalently channels) to detect the activities of primary users. Subsequently, they perform transmission in 

a single channel at variable power and rates depending on the channel sensing decisions and the fading environment. 

A state transition model is constructed to model this cognitive operation. Statistical limitations on the buffer lengths 

y^ • are imposed to take into account the QoS constraints of the cognitive secondary users. Under such QoS constraints 

and limitations on the interference caused to the primary users, the maximum throughput is identified by finding 

^.^ ^ the effective capacity of the cognitive radio channel. Optimal power allocation strategies are obtained and the 

■^ . optimal channel selection criterion is identified. The intricate interplay between effective capacity, interference and 

C^ ' QoS constraints, channel sensing parameters and reliability, fading, and the number of available frequency bands is 

investigated through numerical results. 
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;h ' I. Introduction 



Recent years have witnessed much interest in cognitive radio systems due to their promise as a technology 
that enables systems to utilize the available spectrum much more effectively. This interest has resulted in a 
spur of research activity in the area. In [1], Asghari and Aissa, under constraints on the average interference 
caused at the licensed user over Rayleigh fading channels, studied two adaptation policies at the secondary 
user's transmitter in a cognitive radio system one of which is variable power and the other is variable rate 
and power. They maximized the achievable rates under the above constraints and the bit error rate (BER) 
requirement in MQAM modulation. The authors in ^ derived the fading channel capacity of a secondary 
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user subject to both average and peak received-power constraints at the primary receiver. In addition, they 
obtained optimum power allocation schemes for three different capacity notions, namely, ergodic, outage, 
and minimum-rate. Ghasemi et al. in [3] studied the performance of spectrum- sensing radios under channel 
fading. They showed that due to uncertainty resulting from fading, local signal processing alone may not be 
adequate to meet the performance requirements. Therefore, to remedy this uncertainty they also focused on 
the cooperation among secondary users and the tradeoff between local processing and cooperation in order 
to maximize the spectrum utilization. Furthermore, the authors in [4J focused on the problem of designing 
the sensing duration to maximize the achievable throughput for the secondary network under the constraint 
that the primary users are sufficiently protected. They formulated the sensing-throughput tradeoff problem 
mathematically, and use energy detection sensing scheme to prove that the formulated problem indeed has 
one optimal sensing time which yields the highest throughput for the secondary network. Moreover, Quan et 
al. in [5] introduced a novel wideband spectrum sensing technique, called multiband joint detection, which 
jointly detects the signal energy levels over multiple frequency bands rather than considering one band at a 
time. 

In many wireless systems, it is very important to provide reliable communications while sustaining a 
certain level of quality-of-service (QoS) under time-varying channel conditions. For instance, in wireless 
multimedia transmissions, stringent delay QoS requirements need to be satisfied in order to provide acceptable 
performance levels. In cognitive radio systems, challenges in providing QoS assurances increase due to the 
fact that secondary users should operate under constraints on the interference levels that they produce to 
primary users. For the secondary users, these interference constraints lead to variations in transmit power 
levels and channel accesses. For instance, intermittent access to the channels due to the activity of primary 
users make it difficult for the secondary users to satisfy their own QoS limitations. 

These considerations have led to studies that investigate the cognitive radio performance under QoS 
constraints. Musavian and Aissa in |[6l considered variable-rate, variable-power MQAM modulation em- 
ployed under delay QoS constraints over spectrum- sharing channels. As a performance metric, they used 
the effective capacity to characterize the maximum throughput under QoS constraints. They assumed two 
users sharing the spectrum with one of them having a primary access to the band. The other, known as 
secondary user, is constrained by interference limitations imposed by the primary user. Considering two 
modulation schemes, continuous MQAM and discrete MQAM with restricted constellations, they obtained 
the effective capacity of the secondary user's link, and derived the optimum power allocation scheme that 



maximizes the effective capacity in each case. Additionally, in [|7||, they proposed a QoS constrained power 
and rate allocation scheme for spectrum sharing systems in which the secondary users are allowed to use 
the spectrum under an interference constraint by which a minimum-rate of transmission is guaranteed to the 
primary user for a certain percentage of time. Moreover, applying an average interference power constraint 
which is required to be fulfilled by the secondary user, they obtained the maximum arrival-rate supported 
by a Rayleigh block-fading channel subject to satisfying a given statistical delay QoS constraint. We note 
that in these studies on the performance under QoS limitations, channel sensing is not incorporated into the 
system model. As a result, adaptation of the cognitive transmission according to the presence or absence of 
the primary users is not considered. 

In [[8||, where we also concentrated on cognitive transmission under QoS constraint, we assumed that 
the secondary transmitter sends the data at two different fixed rates and power levels, depending on the 
activity of the primary users, which is determined by channel sensing performed by the secondary users. 
We constructed a state transition model with eight states to model this cognitive transmission channel, and 
determined the effective capacity. On the other hand, we assumed in |[8l that channel sensing is done only 
in one channel, and did not impose explicit interference constraints. 

In this paper, we study the effective capacity of cognitive radio channels where the cognitive radio detects 
the activity of primary users in a multiband environment and then performs the data transmission in one 
of the transmission channels. Both the secondary receiver and the secondary transmitter know the fading 
coefficients of their own channel, and of the channel between the secondary transmitter and the primary 
receiver. The cognitive radio has two power allocation policies depending on the activities of the primary 
users and the sensing decisions. More specifically, the contributions of this paper are the following: 

1) We consider a scenario in which the cognitive system employs multi-channel sensing and uses one 
channel for data transmission thereby decreasing the probability of interference to the primary users. 

2) We identify a state-transition model for cognitive radio transmission in which we compare the trans- 
mission rates with instantaneous channel capacities, and also incorporate the results of channel sensing. 

3) We determine the effective capacity of the cognitive channel under limitations on the average inter- 
ference power experienced by the primary receiver. 

4) We identify the optimal criterion to select the transmission channel out of the available channels and 
obtain the optimal power adaptation policies that maximize the effective capacity. 



5) We analyze the interactions between the effective capacity, QoS constraints, channel sensing duration, 
channel detection threshold, detection and false alarm probabilities through numerical techniques. 

The organization of the rest of the paper is as follows: In Section |Ill we discuss the channel model 
and analyze multi-channel sensing. We describe the channel state transition model in Section |III] under the 
assumption that the secondary users have perfect CSI and send the data at rates equal to the instantaneous 
channel capacity values. In Section |Wl we analyze the received interference power at the primary receiver 
and apply this as a power constraint on the secondary users. In Section |Vl we define the effective capacity 
and find the optimal power distribution and show the criterion to choose the best channel. Numerical results 
are shown in Section |VIl and conclusions are provided in Section IVIII 

II. Cognitive Channel Model and Channel Sensing 

In this paper, we consider a cognitive radio system in which secondary users sense M channels and choose 
one channel for data transmission. We assume that channel sensing and data transmission are conducted 
in frames of duration T seconds. In each frame, N seconds is allocated for channel sensing while data 
transmission occurs in the remaining T — N seconds. Transmission power and rate levels depend on the 
primary users' activities. If all of the channels are detected as busy, transmitter selects one channel with 
a certain criterion, and sets the transmission power and rate to Pk,i{i) and rfci(z), respectively, where 
k E {1, 2, ... , M} is the index of the selected channel and i = 1,2, . . . denotes the time index. Note that if 
Pk,i{i) = 0, transmitter stops sending information when it detects primary users in all channels. If at least 
one channel is sensed to be idle, data transmission is performed with power Pk,2{i) and at rate Tj^ 2(0- If 
multiple channels are detected as idle, then one idle channel is selected again considering a certain criterion. 

The discrete-time channel input-output relation between the secondary transmitter and receiver in the «* 
symbol duration in the /c"^ channel is given by 

yk{i) = hk{i)Xk{i) + nk{i) i = l,2,..., (1) 

if the primary users are absent. On the other hand, if primary users are present in the channel, we have 

yk{i) = hk{i)xk{i) + Sk,p{i) + nk{i) i = l,2,..., (2) 

where Xk{i) and yk{i) denote the complex- valued channel input and output, respectively. In ([T]) and ©, hk{i) 
is the channel fading coefficient between the cognitive transmitter and the receiver. We assume that hk{i) 



has a finite variance, i.e., cr^ < oo, but otherwise has an arbitrary distribution. We define Zk{i) = |/ifc(«)P- 
We consider a block-fading channel model and assume that the fading coefficients stay constant for a block 
of duration T seconds and change from one block to another independently in each channel. In ^, Sk,p{i) 
represents the active primary user's faded signal arriving at the secondary receiver in the A;*^ channel, and 
has a variance a^ (i). nk{i) models the additive thermal noise at the receiver, and is a zero-mean, circularly 
symmetric, complex Gaussian random variable with variance E{|nfc(z)p} = a'^ for all i. We assume that 
the bandwidth of the A;*'* channel is B^. 

In the absence of detailed information on primary users' transmission policies, energy-based detection 
methods are favorable for channel sensing. Knowing that wideband channels exhibit frequency selective 
features, we can divide the band into channels and estimate each received signal through its discrete Fourier 
transform (DFT) [5J. The channel sensing can be formulated as a hypothesis testing problem between the 
noise nk{i) and the signal Sk,p{i) in noise. Noting that there are NBk complex symbols in a duration of 
N seconds in each channel with bandwidth B^, the hypothesis test in channel k can mathematically be 
expressed as follows: 



(3) 



"Hfc.o : yk{i)=nk{i), i = l,...,NBk 
T-Lk,! ■ yk{i) = Sk,p{i) +nk{i), i = l,...,NBk. 
For the above detection problem, the optimal Neyman-Pearson detector is given by [fTOll 

^^ = i^Ei^^»i'^S;:7.. (4) 



NBk 



i=l 



We assume that Sk,p{i) has a circularly symmetric complex Gaussian distribution with zero-mean and variance 
al . Assuming further that {sk,p{i)} are i.i.d., we can immediately conclude that the test statistic Y^ is chi- 
square distributed with 2NBk degrees of freedom. In this case, the probabilities of false alarm and detection 
can be established as follows: 






(5) 



Pk,d = Pr{Yk > 7fe|Hfc,i) = 1 - P I ^^'^' , NBk I (6) 

where P{x, a) denotes the regularized lower gamma function and is defined as P{x, a) = pf^"^ where 7(x, a) 



is the lower incomplete gamma function and T{a) is the Gamma function. In Figure [IJ the probability of 
detection, P^, and the probability of false alarm, Pf, are plotted as a function of the energy detection 
threshold, 7, for different values of channel detection duration. Note that the bandwidth is B = lOkHz 
and the block duration is T = 0.1s. We can see that when the detection threshold is low, Pd and Pf tend 
to be 1, which means that the secondary user, always assuming the existence of an active primary user, 
transmits with power Pi{i) and rate ri(z). On the other hand, when the detection threshold is high, P^ and 
Pf are close to zero, which means that the secondary user, being unable to detect the activity of the primary 
users, always transmits with power P2ii) and rate r2{i), possibly causing significant interference. The main 
purpose is to keep Pd as close to 1 as possible and Pf as close to as possible. Therefore, we have to keep 
the detection threshold in a reasonable interval. Note that the duration of detection is also important since 
increasing the number of channel samples used for sensing improves the quality of channel detection. 

In the hypothesis testing problem in ([3]), another approach is to consider Yk as Gaussian distributed, which 
is accurate if NB^ is large H. In this case, the detection and false alarm probabilities can be expressed in 
terms of Gaussian Q-functions. We would like to note the rest of the analysis in the paper does not depend 
on the specific expressions of the false alarm and detection probabilities. However, numerical results are 
obtained using ^ and ©. 

III. State Transition Model 

In this paper, we assume that both the secondary receiver and transmitter have perfect channel side 
information (CSI), and hence perfectly know the realizations of the fading coefficients {hk{i)}. We further 
assume that the wideband channel is divided into channels, each with bandwidth that is equal to the 
coherence bandwidth Be. Therefore, we henceforth have Bk = Be. With this assumption, we can suppose 
that independent flat fading is experienced in each channel. In order to further simplify the setting, we 
consider a symmetric model in which fading coefficients are identically distributed in different channels. 
Moreover, we assume that the background noise and primary users' signals are also identically distributed 
in different channels and hence their variances a'^ and erf do not depend on k, and the prior probabilities of 
each channel being occupied by the primary users are the same and equal to p. In channel sensing, the same 
energy threshold, 7, is applied in each channel. Finally, in this symmetric model, the transmission power 
and rate policies when the channels are idle or busy are the same for each channel. Due to the consideration 
of a symmetric model, we in the subsequent analysis drop the subscript k in the expressions for the sake 



of brevity. 

First, note that we have the following four possible scenarios considering the correct detections and errors 
in channel sensing: 

Scenario 1: All channels are detected as busy, and channel used for transmission is actually busy. 

Scenario 2: All channels are detected as busy, and channel used for transmission is actually idle. 

Scenario 3: At least one channel is detected as idle, and channel used for transmission is actually 

busy. 

Scenario 4: At least one channel is detected as idle, and channel used for transmission is actually 

idle. 
In each scenario, we have one state, namely either ON or OFF, depending on whether or not the instantaneous 
transmission rate exceeds the instantaneous channel capacity. Considering the interference Sp{i) caused by 
the primary users as additional Gaussian noise, we can express the instantaneous channel capacities in the 
above four scenarios as follows: 

Scenario 1: Ci{i) = BJog^il + SNRi(i)). 

Scenario 2: 02(1) = 5clog2(l + SNR2(i)). 

Scenario 3: C3{i) = 5clog2(l + SNR3(i)). 

Scenario 4: C^^i) = BJog^il + SNR4(i)). 
Above, we have defined 

Pi(i)z(i) , , Pi(i)z(i) , , P2ii)z(i) , , Pi[i)z{i) 

SNRi(2) = } \ . SNR2(^) = \^ \\ SNR3(^) = ^ \ ^ m^^{i) = t,^ - (7) 



Note that z{i) = \h{i)\'^ denotes the fading power. In scenarios 1 and 2, the secondary transmitter detects 
all channels as busy and transmits the information at rate 

ri(^) = 5,log2(l + SNRi(^)). (8) 

On the other hand, in scenarios 3 and 4, at least one channel is sensed as idle and the transmission rate is 

r2« = 5,log2(l + SNR4W), (9) 

since the transmitter, assuming the channel as idle, sets the power level to P2ii) and expects that no 
interference from the primary transmissions will be experienced at the secondary receiver (as seen by the 



absence of a: in the denominator of SNR4). 

In scenarios 1 and 2, transmission rate is less than or equal to the instantaneous channel capacity. Hence, 
reliable transmission at rate ri{i) is attained and channel is in the ON state. Similarly, the channel is in the 
ON state in scenario 4 in which the transmission rate is r2{i). On the other hand, in scenario 3, transmission 
rate exceeds the instantaneous channel capacity (i.e., r2{i) > 6*3(2)) due to miss-detection. In this case, 
reliable communication cannot be established, and the channel is assumed to be in the OFF state. Note that 
the effective transmission rate in this state is zero, and therefore information needs to be retransmitted. We 
assume that this is accomplished through a simple ARQ mechanism. 

For this cognitive transmission model, we initially construct a state transition model. While the ensuing 
discussion describes this model. Figure [2] provides a depiction. As seen in Fig. [21 there are M + 1 ON states 
and 1 OFF state. The single OFF state is the one experienced in scenario 3. The first ON state, which is the 
top leftmost state in Fig. |2l is a combined version of the ON states in scenarios 1 and 2 in both of which 
the transmission rate is ri{i) and the transmission power is Pi{i). Note that all the channels are detected as 
busy in this first ON state. The remaining ON states labeled 2 through (M + 1) can be seen as the expansion 
of the ON state in scenario 4 in which at least one channel is detected as idle and the channel chosen for 
transmission is actually idle. More specifically, the /c* ON state for /c = 2, 3, . . . , M + 1 is the ON state 
in which k — I channels are detected as idle and the channel chosen for transmission is idle. Note that the 
transmission rate is r2{i) and the transmission power is P2{i) in all ON states labeled 2 through (M + 1). 

Next, we characterize the state transition probabilities. State transitions occur every T seconds. We can 

easily see that the probability of staying in the first ON state, in which all channels are detected as busy, is 

expressed as follows: 

Pii = a"" (10) 

where a = pPd + (1 — p) P/ is the probability that channel is detected as busy, and Pd and Pj are the 
probabilities of detection and false alarm, respectively as defined in Q. Recall that p denotes the probability 
that a channel is busy (i.e., there are active primary users in the channel). It is important to note that the 
transition probability in (flOl) is obtained under the assumptions that the primary user activity is independent 
among the channels and also from one block to another. Indeed, under the assumption of independence over 
the blocks, the state transition probabilities do not depend on the originating state and hence we have 

Pii = P21 = ■ ■ ■ = P(M+i)i = P(Af+2)i = a^'^ = Pi (1 1) 



where we have defined pi = pn for alH = 1, 2, . . . , M + 2. Similarly, we can obtain for A; = 2, 3, . . . , M + 1, 



Plk = P2k 



P(M+l)k — P(M+2)k 



p ( (k ~ I) out of M anH ''^^ channel chosen for transmission '\ 

-* Channels are detected as idle '^^^^ is actually idle ) 



(12) 



M 
A;-l 



a^-'+\l-af-'x ^^ ^^^^ ^^^ (13) 

I — a 



■> probability that the channel chosen for 
probabiUty that (fc - 1) out of M channels transmission is actually idle 



are detected as idle 



given that it is detected as idle 



M! 



(M-A; + 1)!(A;-1)! 



a 



M-k+l 



I — a 



>fc-2 



A 
^Pk 



l-p){l-Pf) (14) 
(15) 



Now, we can easily observe that the transition probabilities for the OFF state are 

M+l 
Pl{M+2) = P2(M+2) = • • 



P(M+l){M+2) = P(M+2){M+2) = 1 — 2_^ Plk 

k=l 
M 



(16) 



E 



M! 



^ {M-k)\k\ 



a 



M-k 



^k-l 



I -a) Pil-Pd 



PM+2- 



(17) 



Then, we can easily see that the (M + 2) x (M + 2) state transition probability matrix can be expressed as 



R 



Pi,i 



Pi, M+2 



Pi 



PM+2 



PM+2,1 ■ ■ PM+2,M+2 Pi ■ ■ Pm+2 

Note that R has a rank of 1. Note also that in each frame duration of T seconds, ri(k)(T — N) bits are 
transmitted and received in state 1, and r2{k)(T — A^) bits are transmitted and received in states 2 through 
M + l, while the transmitted number of bits is assumed to be zero in state M + 2. 



IV. Interference Power Constraints 

In this section, we consider interference power constraints to limit the transmission powers of the secondary 
users and provide protection to primary users. In particular, we assume that the transmission power of the 
secondary users is constrained in such a way that the average interference power on the primary receiver is 
limited. 



Note that interference to the primary users is caused in scenarios 1 and 3. In scenario 1, the channel 
is busy, and the secondary user, detecting the channel as busy, transmits at power level Pi. Consequently, 
the instantaneous interference power experienced by the primary user is PiZgp where Zsp = \hsp{i)\'^ is the 
magnitude-square of the fading coefficient of the channel between the secondary transmitter and the primary 
user. Note also that the probability of being in scenario 1 (i.e., the probability of detecting all channels busy 
and having the chosen transmission channel as actually busy) is a^^^^pP^, as can be easily seen through 
an analysis similar to that in (fT3l) . 

In scenario 3, the secondary user, detecting the channel as idle, transmits at power P2 although the 
channel is actually is busy. In this case, the instantaneous interference power is P2Zsp- Since we consider 
power adaption, transmission power levels Pi and P2 in general vary with Zgp and also with z, which is 
the power of the fading coefficient between the secondary transmitter and secondary receiver in the chosen 
transmission channel. Hence, in both scenarios, the instantaneous interference power levels depend on both 
Zsp and z whose distributions depend on the criterion with which the transmission channel is chosen and 
the number of available channels from which the selection is performed. For this reason, it is necessary in 
scenario 3 to separately consider the individual cases with different number of idle-detected channels. We 
have M such cases. For instance, in the A;*'^ case for k = 1,2,..., M, we have k channels detected as idle 
and the channel chosen out of these k channels is actually busy. The probability of the A;*^ case can be 
easily found to be j^^^ a^'-^ (1 - af'^ p(l - P,). 

Following the above discussion, we can now express the average interference constraints as follows: 

^ M' 

«*^^V^ E{PiZsp\ +I]7^^^^r^y|^«*'^'(l-«)'~V(l-^d) E^{P2Zsp} <h,g (18) 

probability of average interference ~ ^ v '' average interference 

scenario 1 in scenario 1 probability of the fc*'' case of scenario 3 in the fc"' case 

of scenario 3 

Note from above that lavg is the constraint on the interference averaged over the distributions of z and 
Zgp (through the expectations), and also averaged over the probabilities of different scenarios and cases. It 
is important to note that the term Ek {P2Zsp}, as discussed above, depends in general on the number of 
idle-detected channels, k. This dependence is indicated through the subscript k. 

In a system with more strict requirements on the interference, the following individual interference 
constraints can be imposed: 

E{PiZsp}<Io and Et{P2Zsp]<h for fc = 1, 2, . . . , M. (19) 

10 



If, for instance, Iq = Ii = I2 = ■ ■ ■ = Im, then interference averaged over fading is limited by the same 
constraint regardless of which scenario is being realized. As considered in Q, by appropriately choosing the 
values of Jo and Ik in (fT9l ), we can provide primary users a minimum rate guarantee for a certain percentage 
of the time in a Rayleigh fading environment through the following outage constraints: 



Pr <! log2 I 1 + — -%|^^^L^ I < ^™™ r^ < ^2T- forfc = l,2,...,M. (21) 






Lp 



P°"* and Pg^fc* '^^^ be seen as the outage constraints in scenario 1 and in the k^'^ case of scenario 3, 
respectively. In the above formulations, Rmin is the required minimum transmission rate to be provided 
to the primary users with outage probabilities P°"* and P2^, and Zp(i) = \hp(i)\'^ where hp is the fading 
coefficient of the channel between the primary transmitter and primary receiver, a^ is the variance of 
the zero-mean, circularly symmetric, complex Gaussian thermal noise at the primary receiver. Ppri is the 
transmission power of the primary transmitter. Under the assumption that Zp is an exponential random 
variable (i.e., we have a Rayleigh fading channel between the primary transmitter and receiver), the outage 
probability in (|20|) can be expressed as follows: 



^^ {'"^'^ {' + p.( wl; b) ^ «"■"•} - '^ {- ^ ^ (^-("^-W + <^0} <^^' 



-^^=^^i-i(PiW..p«+<B,) 



--l{E{P,{i)z,pit)}+al^B,) 



<l-e ^p" v-i-vv^^p^vj."np^c; ^24) 

where (1231) is obtained by performing integration with respect to the probability density function (pdf) of 
Zp in the evaluation of the probability expression in (|22l) . As a result, the expectation in (l23l) is with respect 
to the remaining random components Pi and Zgp. Finally, the inequality in (|24)) follows from the concavity 
of the function 1 — e^^ and Jensen's inequality. From (|24l) . we can immediately see that if we impose 



out\ 



E {P,z,p} < $1 = -^4,-—L^ - <P„ (25) 



2, TTiin —1 
p '. 

then the constraint in (|20l) will be satisfied. A similar discussion follows for (|2TI) as well. 

In the subsequent parts of the paper, we assume that an average interference power constraint in the form 

11 



given in (1181 ) is imposed. 



V. Effective Capacity 



In this section, we identify the maximum throughput that the cognitive radio channel with the aforemen- 
tioned state-transition model can sustain under interference power constraints and statistical OoS limitations 

n 

imposed in the form of buffer or delay violation probabilities!!. Wu and Negi in [|TT| defined the effective 
capacity as the maximum constant arrival rate that can be supported by a given channel service process 
while also satisfying a statistical QoS requirement specified by the QoS exponent 9. If we define Q as the 
stationary queue length, then 9 is defined as the decay rate of the tail distribution of the queue length Q: 

li,„ !2i^WM = _,^ (26) 

g-^oo q 

Hence, we have the following approximation for the buffer violation probability for large Qmax'- P{Q > 
Qmax) ~ e"^''"""'. Therefore, larger 9 corresponds to more strict QoS constraints, while the smaller 9 implies 
looser constraints. In certain settings, constraints on the queue length can be linked to limitations on the delay 
and hence delay-QoS constraints. It is shown in [12J that P{D > dmax} < c^/P{Q > Qmax} for constant 
arrival rates, where D denotes the steady-state delay experienced in the buffer. In the above formulation, c 
is a positive constant, Qmax = cidmax and a is the source arrival rate. Therefore, effective capacity provides 
the maximum arrival rate when the system is subject to statistical queue length or delay constraints in the 
forms of P{Q > gmax) < e^^'^™"^ or P{D > t/max} < ce^^""^™"^/^, respectively. Since the average arrival 
rate is equal to the average departure rate when the queue is in steady-state [|T3l . effective capacity can also 
be seen as the maximum throughput in the presence of such constraints. 
The effective capacity for a given QoS exponent 9 is given by 

-limllog,E{e-^^W} = -^ (27) 

where S{t) = J2k=i ^(^) ^^ ^^^ time-accumulated service process, and {r(k), A; = 1, 2, . . . } is defined as the 
discrete-time, stationary and ergodic stochastic service process. Note that A{9) is the asymptotic log-moment 

'Note that interference constraints are imposed to provide a certain level of quality-of-service to the primary users, while buffer or delay 
constraints are used to statistically guarantee a quality-of-service level to the transmissions of the secondary users. Hence, the formulation in 
the paper effectively considers service guarantees for both the primary and secondary users. On the other hand, QoS constraints throughout the 
paper refer to buffer/delay constraints to avoid confusion. 
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generating function of S{t), and is given by 

A(^) = lim -logETe^'^^*)] . (28) 

The service rate according to the model described in Section [ni] is r(k) = ri(k){T — N) if the cognitive 
system is in state 1 at time k. Similarly, the service rate is r{k) = r2{k){T — N) in the states between 2 
and M + 1. In the OFF state, instantaneous transmission rate exceeds the instantaneous channel capacity 
and reliable communication can not be achieved. Therefore, the service rate in this state is effectively zero. 

In the next result, we provide the effective capacity for the cognitive radio channel and state transition 
model described in the previous section. 

Theorem 1: For the cognitive radio channel with the state transition model given in Section [nil the 
normalized effective capacity (in bits/s/Hz) under the average interference power constraint (flSl) is given by 

R^iSNR, 9) = --i- max log, (p,E {e-(^-^)^- } + f^p.+^E, {e-^T-N)er. | ^ p\ 

_:::J-avg 

(29) 

Above, pk for A; = l,2,...,M + 2 denote the state transition probabilities defined in (fTTI) . (fT5l) . and (fTTI) in 
Section Unl Note also that the maximization is with respect to the power adaptation policies Pi and P2. 

Remark: In the effective capacity expression (|29l) . the expectation E {PiZgp} in the constraint and E |e~'^^~^^^''^ } 
are with respect to the joint distribution of {z, Zsp) of the channel selected for transmission when all channels 
are detected busy. The expectations E^ {P2Zsp} and E^ ^e~^'^~^^^^'^ } are with respect to the joint distribution 
of {z, Zsp) of the channel selected for transmission when k channels are detected as idle. 

Proof of TheoremU} In [|9i Chap. 7, Example 7.2.7], it is shown for Markov modulated processes that 

!^ = ^\og^sp{mR) (30) 

where sp{(p{9)R) is the spectral radius (i.e., the maximum of the absolute values of the eigenvalues) of the 
matrix 0(6')_R, R is the transition matrix of the underlying Markov process, and (f)(9) = diag(0i(^), . . . , 0m+2(^)) 
is a diagonal matrix whose components are the moment generating functions of the processes in given states. 
The rates supported by the cognitive radio channel with the state transition model described in the previous 
section can be seen as a Markov modulated process and hence the setup considered in |91 can be immediately 
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applied to our setting. Since the processes in the states are time-varying transmission rates, we can easily 
find that 0(e) = diag {E {e^^-^^^^^ } , E^ {e(T-NWr2 },..., Em {e^^-^)^'-^ } , l}. Then, we have 



(t){9)R 



h{9)pi 



(f)l{0)PM 



+2 



4'M+2{d)Pl ■ ■ (f)M+2{d)PM+2 

Since (j){9)R is a matrix with unit rank, we can readily find that 



sp{4>{e)R) = trace (^0(e)i?j = 0i(e)pi + 02(^)P2 + ■ ■ ■ + 0m+i(^W+i + 0m+2(^K/+2 (31) 

= p^E {e(^-^)^- } + p,E, {e(^-^)^- } + ■■■+ pm-^iE^ {e^^"^)^- } + pm+2. (32) 



Then, combining (1321) with (1301) and (|27l) . normalizing the expression with TBc in order to have the effective 
capacity in the units of bits/s/Hz, and considering the maximization over power adaptation policies, we reach 
to the effective capacity formula given in (|29l) . D 

We would like to note that the effective capacity expression in ( |29l) is obtained for a given sensing 
duration N, detection threshold 7, and QoS exponent 6. In the next section, we investigate the impact of 
these parameters on the effective capacity through numerical analysis. Before the numerical analysis, we 
first identify below the optimal power adaptation policies that the secondary users should employ. 

Theorem 2: The optimal power adaptations for the secondary users under the constraint given in (fTSi) are 



Ml 



Pi 



Zspfi\\ 



1 

c+1 



and 



M2 



P9 



Zap 1^2^ 



1 
c+1 



otherwise 



^>/32A 
otherwise 



(33) 



where /^i = i?c(cr^ + c^l 



/^2 



alB,, c = B,{T - N)9/ log, 2, /3i = ^ and (32 



Pil-Pd)f^2 



(34) 



A is a 



c(l-p)(l-P^)- 

parameter whose value can be found numerically by satisfying the constraint (fTSi) with equality. 

Proof: Since logarithm is a monotonic function, the optimal power adaptation policies can also be obtained 
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from the following minimization problem: 

M 

min p,E {e-(^-^)^- } + Vp.+iE, {e-(^-^)^- } (35) 



It is clear that the objective function in (1351) is strictly convex and the constraint function in (1181) is linear with 
respect to Pi and P2 Q Then, forming the Lagrangian function and setting the derivatives of the Lagrangian 
with respect to Pi and P2 equal to zero, we obtain: 

-c-l" 

a^'f{z,Zsp) = Q (36) 



/^2 V /^2 / 



-c-l 



M! 



5^ ^^^-'^(l - ^r'^M^^i^M^^ ^s,) = (37) 

where A is the Lagrange multiplier. Above, f{z, Zgp) denotes the joint distribution of {z, Zgp) of the channel 
selected for transmission when all channels are detected busy. Hence, in this case, the transmission channel is 
chosen among M channels. Similarly, fk{z, Zsp) denotes the joint distribution when k channels are detected 
idle, and the transmission channel is selected out of these k channels. Defining /3i = i^i£^ and /32 = 
c(i )(i p) ' ^^^ solving (l36l) and (l37l) . we obtain the optimal power policies given in (|33l) and (|34l) . D 
Now, using the optimal transmission policies given in (|33l) and (|34l) . we can express the effective capacity 
as follows: 

Above, the subscripts /3iA and /32A in the expectations denote that the lower limits of the integrals are equal 
these values and not to zero. For instance, E^^x { ( ^ ^b \ 1 '^ \ ~ //Ta ( W\ ) " /^(•^) '^•^• 

Until now, we have not specified the criterion with which the transmission channel is selected from a 
set of available channels. In (l38l) . we can easily observe that the effective capacity depends only on the 

c 

channel power ratio ^, and is increasing with increasing -j- due to the fact that the terms f ^ ^ ^ j 

— c 

and [ g ,^ j '' are monotonically decreasing functions of j-. Therefore, the criterion for choosing the 
transmission band among multiple busy bands unless there is no idle band detected, or among multiple idle 

^Strict convexity follows from the strict concavity of r\ and r2 in ([8} and ([9} with respect to Pi and P2 respectively, strict convexity of the 
exponential function, and the fact that the nonnegative weighted sum of strictly convex functions is strictly convex 1141 Section 3.2.1]. 
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bands if there are idle bands detected should be based on this ratio of the channel gains. Clearly, the strategy 
that maximizes the effective capacity is to choose the channel (or equivalently the frequency band) with the 
highest ratio of ^. This is also intuitively appealing as we want to maximize z to improve the secondary 
transmission and at the same time minimize Zsp to diminish the interference caused to the primary users. 
Maximizing -^ provides us the right balance in the channel selection. 

We define x = maxjgn 2 ... m} -^^ where -^^ is the ratio of the gains in the i^^ channel. Assuming that 
these ratios are independent and identically distributed in different channels, we can express the pdf of x as 

(39) 



Ux) = Mf^{x) 



Zsp 



Af-1 



where /^ and F^_ are the pdf and cumulative distribution function (cdf), respectively, of ^, the gain 



c + l 



ratio in one channel. Now, the expectation E^^a \ [ z ^0 x ) f"' which arises under the assumption that 

all channels are detected busy and the transmission channel is selected among these M channels, can be 
evaluated with respect to the distribution in (|39l ). 

Similarly, we define Xk = maxjgn 2 k} -^^ for k = 1, . . . , M. The pdf of Xk can be expressed as follows: 



fx.ix) = kfMx) 



F^(x) 



fc-i 



k = l,2,...,M. 



(40) 



The expectation E^^p^x 



Zsp&A 



can be evaluated using the distribution in ( |40l ). Finally, after some 



calculations, we can write the effective capacity in integral form as 



Re {SNR, I 



A/-1 






:i-p)il-Pj)M f^ix) a + il-a)F^ix) 



X 

M-1 



c + l 



dx 
/92A 



X 



c+l 



dx+PM+2>- (41) 



VI. Numerical Results 

In this section, we present numerical results for the effective capacity as a function of the channel sensing 
reliability (i.e., detection and false alarm probabilities) and the average interference constraints. Throughout 
the numerical results, we assume that QoS parameter is 9 = 0.1, block duration is T = Is, channel sensing 
duration is A^ = 0.1s, and the prior probability of each channel being busy is p = 0.1. 

Before the numerical analysis, we first provide expressions for the probabilities of operating in each one 
of the four scenarios described in Section Unl These probabilities are also important metrics in analyzing 
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the performance. We have 
Pjsecondary system is in scenario 1} = Ps^ = a^'^^pPd, 

Pjsecondary system is in scenario 3} = P53 = Y] a*^"^(l — aY — 



Pjsecondary system is in scenario 2} = Pg^ = a*^ ^(1 ~ p)Pf 

M 



fc=l 




> probability that tlie channel chosen 
probability that at least one channel f™' transmission is actually busy 

is detected as idle g'™" *at it is detected as idle 



(42) 



^_(l-a^Xl-P, 
P{ secondary system is in scenario 4} = Ps^ 



1 — a 
(l-a*0(l-p)(l-P;) 



1 — a 

In Figure |3l we plot these probabilities as a function of the detection probability Pd for two cases in 
which the number of channels is M = 1 and M = 10, respectively. As expected, we observe that P^^ and 
P52 decrease with increasing M. We also see that P53 and P^^ are assuming small values when P^ is very 
close to 1. Note from Fig. [T]that as Pd approaches 1, the false alarm probability P/ increases as well. 

A. Rayleigh Fading 

The analysis in the preceding sections apply for arbitrary joint distributions of z and Zsp under the mild 
assumption that the they have finite means (i.e., fading has finite average power). In this subsection, we 
consider a Rayleigh fading scenario in which the power gains z and Zsp are exponentially distributed. We 
assume that z and Zsp are mutually independent and each has unit-mean. Then, the pdf and cdf of -^ can 
be expressed as follows: 

f^(x) = x>0 and P^fx) = — ^ x > 0. (43) 

^"p ^ (x + 1)2 - -"P ^ ^ X + 1 

In Fig. m we plot the effective capacity vs. probability of detection, P^, for different number of channels 
when the average interference power constraint normalized by the noise power is Iavg{dB) = 10 logio ( 2°"o ) 
OdB, where a.^ is the noise variance at the primary user. We observe that with increasing Pa, the effective 
capacity is increasing due to the fact more reliable detection of the activity primary users leads to fewer 
miss-detections and hence the probability of scenario 3 or equivalently the probability of being in state 
M + 2, in which the transmission rate is effectively zero, diminishes. We also interestingly see that the 
highest effective capacity is attained when M = 1. Hence, secondary users seem to not benefit from the 
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availability of multiple channels. This is especially pronounced for high values of P^. Although several 
factors and parameters are in play in determining the value of the effective capacity, one explanation for 
this observation is that the probabilities of scenarios 1 and 2, in which the secondary users transmit with 
power Pi, decrease with increasing M, while the probabilities of scenarios 3 and 4 increase as seen in (|42)) . 
Note that in scenario 3, no reliable communication is possible and transmission rate is effectively zero. In 
Fig. |5l we display similar results when lavg = — lOdB. Hence, secondary users operate under more stringent 
interference constraints. In this case, we note that M = 2 gives the highest throughput while the performance 
with M = 1 is strictly suboptimal. 

In Fig.[6l we show the effective capacities as a function lavg (dB) for different values of M when Pd = 0.9 
and Pf = 0.2. Confirming our previous observation, we notice that as the interference constraint gets more 
strict and hence lavg becomes smaller, a higher value of M is needed to maximize the effective capacity. 
For instance, M = 10 channels are needed when lavg < — 30dB. On the other hand, for approximately 
lavg > — 6dB, having M = 1 gives the highest throughput. 

Above, we have remarked that increasing the number of available channels from which the transmission 
channel is selected provides no benefit or can even degrade the performance of secondary users under certain 
conditions. On the other hand, it is important to note that increasing M always brings a benefit to the primary 
users in the form of decreased probability of interference. In order to quantify this type of gain, we consider 
below the probability that the channel selected for transmission is actually busy and hence the primary user 
in this channel experiences interference: 

„ T-. /channel selected\ ^ /channel selected , ,, ^v,ar,„orc n™\ x^ /channel selected , , r ^ u r\ ,. .-. 

p.. = P for transmission = P for transmission and £tpJ?P^" i hfs„ + P for transmission and ^^ '^^f °"5 '^''^f ' (44) 

^'^^ V is actually busy J \k actually busy detected as busy^ "T -^ ^ .^ ^^^^^^jj^ ^^^^ aiiu ^^^ detected as idle J \^^J 

= Ps, + Ps, (45) 

l-a'^-P, + P,a''-' 

= P -. ■ (46) 

I — a 

Note that Pint depends on P^ and also Pf through a = pPd + (1 — p)Pf. It can be easily seen that this 
interference probability Pint decreases with increasing M when Pa > Pf. As M goes to infinity, we have 
limM^oo -Pint = P^jz^- Indeed, in this asymptotic regime. Pint becomes zero with perfect detection (i.e., 
with Pd = I). Note that secondary users transmit (if Pi > 0) even when all channels are detected as busy. 
As M — 7- oo, the probability of such an event vanishes. Also, having P^ = 1 enables the secondary users to 
avoid scenario 3. Hence, interference is not caused to the primary users. 
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In Fig. |71 we plot Pint vs. the detection probability for different values of M. We also display how the 
false alarm probability evolves as P^ varies from to 1. It can be easily seen that while Pint = p when 
M = 1, a smaller Pint is achieved for higher values of M unless P^ = 1. On the other hand, as also 
discussed above, we immediately note that Pint monotonically decreases to as Pd increases to 1 when M 
is unbounded (i.e., M — )■ oo). 

B. Nakagami Fading 

Nakagami fading occurs when multipath scattering with relatively large delay-time spreads occurs. There- 
fore, Nakagami distribution matches some empirical data better than many other distributions do. With this 
motivation, we also consider Nakagami fading in our numerical results. The pdf of the Nakagami-m random 
variable y = \h\ is given by fy{y) = f^ ( ^ ) y'^"^~^e ^"y where m is the number of degrees of freedom. 
If both Zgp and z have the same number of degrees of freedom, we can express the pdf of a; = ^ as follows: 



^sp 



Note also that Rayleigh fading is a special case of Nakagami fading when m = 1. In our experiments, we 
consider the case in which m = 3. Now, we can express the cdf of x for m = 3 as 

[X + 1)^ [x + ly^ [x + 1)^ 

In Fig. [HI we plot effective capacity vs. lavg (dB) for different values of M when P^ = 0.9 and Pj = 0.2. 
Here, we again observe results similar to those in Fig. [6l We obtain higher throughput by sensing more than 
one channel in the presence of strict interference constraints on cognitive radios. 

VII. Conclusion 

In this paper, we have studied the performance of cognitive transmission under QoS constraints and 
interference limitations. We have considered a scenario in which secondary users sense multiple channels 
and then select a single channel for transmission with rate and power that depend on both sensing decisions 
and fading. We have constructed a state transition model for this cognitive operation. We have meticulously 
identified possible scenarios and states in which the secondary users operate. These states depend on sensing 
decisions, true nature of the channels' being busy or idle, and transmission rates being smaller or greater 
than the instantaneous channel capacity values. We have formulated and imposed an average interference 
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constraint on the secondary users. Under such interference constraints and also statistical QoS limitations 
in the form of buffer constraints, we have obtained the maximum throughput through the effective capacity 
formulation. Therefore, we have effectively analyzed the performance in a practically appealing setting in 
which both the primary and secondary users are provided with certain service guarantees. We have determined 
the optimal power adaptation strategies and the optimal channel selection criterion in the sense of maximizing 
the effective capacity. We have had several interesting observations through our numerical results. We have 
shown that improving the reliability of channel sensing expectedly increases the throughput. We have noted 
that sensing multiple channels is beneficial only under relatively strict interference constraints. At the same 
time, we have remarked that sensing multiple channels can decrease the chances of a primary user being 
interfered. 
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Fig. 2. State transition model for the cognitive radio channel. The numbered label for each state is given on the lower-right corner of the box 
representing the state. 
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